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ABSTRACT 

Based on the radiation hydrodynamical model for the black hole (BH) growth, incorporated with the 
chemical evolution of the early-type host galaxy, we construct the coevolution model of a QSO BH and 
the host galaxy. As a result, it is found that after a galactic wind epoch, the luminosity is shifted from 
the host-dominant phase to the AGN-dominant phase (QSO phase) in the timescale of a few 10^ years. 
The former phase corresponds to the early stage of growing BH, and can be regarded as a "proto-QSO" 
phase. It has observable characteristic properties as follows: (1) The width of broad emission line is 
narrower than that of ordinary QSOs, and it is typically less than 1500km/s. (2) The BH-to-bulge mass 
ratio, MBH/-^^buige, is in the range of 10~^'^ — 10~'^'^. (3) Host galaxies are bluer compared to QSO 
hosts, by about 0.5 magnitude in the colors of {B — V) at the rest bands and {V — K) at the observed 
bands, with assuming galaxy formation redshifts of Zf = 3 — 5. (4) The metallicity of gas in galactic 
nuclei is ~ 8Zq, and that of stars weighted by the host luminosity is ^ 3Zq. (5) The central massive 
BH (~ 1O^M0) is surrounded by a massive dusty disk (> IO^AIq), which may obscure the nucleus in 
the edge-on view and make a type 2 nucleus. By comparing these predictions with recent observations, 
radio galaxies are a possible candidate for proto-QSOs. Also, it is anticipated that the proto-QSO phase 
is preceded by an optically thick phase, which may correspond to ULIRGs. In this phase, A/bh /-^^buigo 
is predicted to be much less than 10^'^ and grow with metallicity. Moreover, as precursors of ULIRGs, 
optically-thin star-forming galaxies are predicted. These may be in the assembly phase of Lyman break 
galaxies (LBGs) or Lya emitters. 

Subject headings: black hole physics — galaxies: active — 

galaxies : evolution — galaxies : formation — galaxies : nuclei — galaxies : 
starburst — radiation drag 



1. INTRODUCTION 

Recent X-ray and optical observations suggest the pos- 
sibility that active galactic nuclei (AGNs) could be divided 
into two subclasses according to the rate of black hole (BH) 
growth; one is a rapidly growing phase and the other is 
a slow growing phase (Pounds et al. 1995; Boiler et al. 
1996; Mineshige et al. 2000; Mathur et al. 2001: Wandel 
2002). This possibility has been pointed out primarily for 
the Seyfert 1 galaxies (Syls), which are divided into two 
subclasses according to the width of broad emission line, 
%LR- Syls with Vblr less than 2000km/s are called nar- 
row line Syls (NLSyls), whereas those with broader line 
width are called broad hne Syls (BLSyls). NLSyls ex- 
hibit two distinctive X-ray properties, that is, rapid X-ray 
variability and strong soft X-ray excess. These properties 
can be explained in terms of the optically thick ADAF 
(advection-dominated accretion flow) onto a smaller BH, 
which is realized by higher accretion rate compared to the 
Eddington hmit (Pounds et al. 1995; Boiler et al. 1996; 
Mineshige et al. 2000). Also, it is pointed out that the BH- 
to-bulge mass ratio is noticeably smaller than that in el- 
hptical galaxies, MeH/Mbuige < 10"^ (Mathur et al. 2001; 
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Wandel 2002). All of these suggest that NLSyls are in the 
rapidly growing phase of BH, in contrast to BLSyls which 
are explained by conventional mild accretion onto a large 
BH. Additionally, Kawaguchi & Aoki (2000) argue that 
NLSyls have a high star formation rate (SFR). According 
to these observations, it has been suggested that NLSyls 
may be Syls in the early stage of their evolution (Mathur 
2000). On the analogy of NLSyls, QSOs are also expected 
to have rapidly growing phase of QSO BHs. But, it has 
not been elucidated so far what objects correspond to the 
early phase of QSOs. 

On the other hand, recent high-resolution observations 
of galactic centers have revealed that the estimated mass 
of a central "massive dark object" (MDO), which is the 
nomenclature for a supermassive BH candidate, does cor- 
relate with the mass of a galactic bulge; the mass ratio of 
the BH to the bulge is 0.001-0.006 as a median value (Ko- 
rmendy & Richstone 1995; Richstone et al. 1998; Magor- 
rian et al. 1998; Loar 1998; Gebhardt et al. 2000; Fer- 
rarese & Merritt 2000; Merritt & Ferrarese 2001 McLure 
& Dunlop 2001; McLure & Dunlop 2002; Wandel 2002). 
(It is noted that the bulge means a whole galaxy for an 
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elliptical galaxy.) In addition, it has been found that QSO 
host galaxies are mostly luminous and well-evolved early- 
type galaxies (McLeod & Rieke 1995; Bahcall et al. 1997; 
Hooper, Impey & Foltz 1997; McLoed, Rieke fc Storrie- 
Lombardi 1999; Brotherton et al. 1999; Kirhakos et al. 
1999; McLure et al. 1999; McLure, Dunlop & Kukula 

2000) . These findings, combined with the BH-to-bulge 
relations, suggest that the formation of a supermassive 
BH, an elliptical galaxy, and a QSO is physically related 
to each other. But, the link between the formation of a 
supermassive BH and the evolution of a host galaxy is 
an open question. Also, the physical relationship among 
QSOs, ultraluminous infrared galaxies (ULIRGs), and ra- 
dio galaxies has been an issue of long standing. 

Some theoretical models of BH growth models have 
been considered to explain the BH-to-bulge correlations 
(Silk & Rees 1998; Ostriker 2000; Adams, Graff, & Rich- 
stone 2001). But, little has been elucidated regarding the 
physics on the angular momentum transfer, which is requi- 
site for BH formation. Recently, as a potential mechanism 
to remove angular momentum, Umemura (2001) has con- 
sidered the effects of radiation drag, which is equivalent 
to a well-known Poynting- Robertson effect. The exact ex- 
pressions for the radiation drag are found in the literature 
(Umemura, Fukue, & Mineshige 1997; Fukue, Umemura, 
& Mineshige 1997). In an optically thick regime, the effi- 
ciency of radiation drag is saturated due to the conserva- 
tion of the photon number (Tsuribe, & Umemura 1997). 
Thus, the angular momentum loss rate by the radiation 
drag is given by rfln J/dt ~ — (L*/c^Mg), where J, Z/*, and 
Mg are the total angular momentum of gaseous compo- 
nent, the total luminosity of the bulge, and the total mass 
of gas, respectively. Then, the maximal rate of mass ac- 
cretion is given by M = — Mgdin J/dt = i*/c^ (Umemura 

2001) . Thus, the total accreted mass on to the MDO, 
-^MDO, is estimated by 

Mmdo^ / -^dt. 
Jo c 

In practice, the interstellar medium (ISM) is observed to 

be highly inhomogeneous in an active star-forming galax- 
ies (Sanders et al. 1988; Gordon, Calzetti & Witt 1997). 
Kawakatu & Umemura (2002) have shown that the inho- 
mogeneity of interstellar medium helps the radiation drag 
to sustain the maximal efficiency. 

Thus, the final mass of MDO is proportional to the to- 
tal radiation energy from bulge stars, and the resultant 
BH-to-bulge mass ratio is basically determined by the en- 
ergy conversion efficiency of the nuclear fusion from hy- 
drogen to helium, i.e. 0.007 (Umemura 2001). So far, 
the realistic chemical evolution of the host galaxy has not 
been incorporated, but a simple evolutionary model was 
assumed. As for the relation between a QSO BH and 
the host galaxy, some phenomenological models have been 
proposed (Haehnelt & Rees 1993; Haiman & Loeb 1998; 
Kauffmann & Haehnelt 2000; Monaco, Salucci, & Danese 
2000; Granato et al. 2001; Hosokawa et al. 2001), but 
little on the physics has been known. Hence, in order 
to reveal the formation and evolution of QSOs and clar- 
ify what objects correspond to the early phase of QSOs, 
it is important to investigate the physics on the rapidly 
growing phase of QSO BHs. Here, based on the radiation 
drag model with incorporating the realistic chemical evo- 



lution, we construct a physical model for the coevolution 
of a QSO BH and the early-type host galaxy. The pur- 
pose of this paper is to elucidate the physical relationship 
between a BH growth and the evolution of host galaxy, 
and define a proto-QSO phase as an early stage of QSO 
evolution. Then, we predict the observable properties of 
proto-QSOs. Also, we address a unified picture for the 
evolution of an elliptical galaxy nucleus. 

The paper is organized as follows. In section 2, we build 
up a theoretical model for the coevolution of a QSO BH 
and the early- type host galaxy. In section 3, we investi- 
gate the time-dependent relation between a QSO BH and 
the early-type host galaxy, and analyze the physical states 
of proto-QSOs which correspond to the rapidly growing 
phase of a QSO BH. In Section 4, we propose a unified 
picture for the evolution of elliptical galaxy nucleus. Sec- 
tion 5 is devoted to the conclusions. 



2. COEVOLUTION MODEL 

First, we model the BH growth based on the radi- 
ation drag-driven mass accretion. Here, we suppose a 
two-component system that consists of a spheroidal stel- 
lar bulge and inhomogeneous optically-thick interstellar 
medium within it. The radiation drag efficiency increases 
with the optical depth t in proportion to (1 — e""^) 
(Umemura 2001). Thus, the mass of an MDO, Muuo, 
which is the total mass of dusty ISM assembled to the 
central massive object, is given by 

Muuoit) = 7?drag f f b^^^l^ (l - e--(*)) dudt, 

(1) 

where t,^ is defined as the optical depth of the bulge mea- 
sured from the center. ?7drag is found to be maximally 
0.34 in the optically thick regime (Kawakatu & Umemura 
2002). Here, we estimate the evolution of for U-, B-, 
V- and K-hand by using an evolutionary spectral synthe- 
sis code 'PEGASE' (Fioc & Rocca-Volmerange 1997). In 
Figure 1, the evolution of Tj, for all bands is shown. We 
define a time tthin before which the optical depth is less 
than unity in the ?7-band, which is most intensive in the 
early evolutionary phase. Before tthin, the radiation drag 
efficiency is low and the growth of MDO is quite slow. 
Furthermore, once a galactic wind occurs, the bulge be- 
comes ISM-deficient and thus optically thin, so that the 
mass accretion via radiation drag is terminated. In this 
paper, after a wind epoch t^, of several 10*yr, is as- 
sumed to drop abruptly to r,^ ^ 1. Hence, the final mass 
of an MDO is given by 

MMDO = %rag r ^^^^rft, (2) 

J tthln ^ 

where Lbuige is the bolometric luminosity of the bulge. It 
is found that escaping photons from the bulge in the opti- 
cally thin phase ofO <t < fthin is less than 5% of the total. 
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Fk;. 1. The optical depth (r^) of bulge measured from the 
center as a function of time. The abscissa is time in units of yr. 
The ordinate is the optical depth for U—,B — , V—, and K— band. 
A time tthin is defined so that the optical depth is less than unity 
in the (/-band. The optically thin phase (rjj < 1) is the gray area. 
The system is optically thick until the galactic wind timcscale (tw), 
at which Ti, is assumed to drop abruptly down to Tv <SC 1. 

In this model, we should distinguish the BH mass from 
the mass of an MDO, although the mass of an MDO is 
often regarded as BH mass from an observational point 
of view. Supposing the mass accretion driven by the vis- 
cosity on to the BH horizon is determined by an order of 
Eddington rate, the BH mass grows according to 

Mbh = Moe''*/*^■^^ (3) 

where v is the ratio of BH accretion rate to the Edding- 
ton rate, — Mbk/ M-^dd, which is about 0.1 for QSOs 
(McLeoed, Rieke,& Storrie-Lombaardi 1999), and tEdd is 
the Eddington timescale, tEdd = VbhMbhc^ / Lsdd, where 
r]BH is the energy conversion efficiency and L^dd is the 
Eddington luminosity. Recently, it is shown by a full 
general rclativistic calculation (Shibata & Shapiro 2002) 
that the collapse of a rotating supermassive star (SMS) 
results in the formation of a rotating BH with Kerr pa- 
rameter of 0.75. Hence, 77drag is assumed to be 0.42, 
which is the efficiency of a extreme Kerr BH. Then, we 
have tEdd = 1.9 X lO^yr. Mq is the mass of a seed BH, 
which could be an early formed massive BH with ~ lO^M© 
(Umemura, Loeb, & Turner 1993) or a massive BH with 
^ 10"'' Mq formed by the collapse of a rotating supermas- 
sive star (SMS) (Baumgarte & Shapiro 1999; Saijo et al. 
2002; Shibata 2002), which may result from by the viscous 
runaway collapse of the MDO (Tsuribe 1999; Umemura 
2002). 

Next, we construct the model for the chemical evolution 
of host galaxy. To treat the realistic chemical evolution, 
we use an evolutionary spectral synthesis code 'PEGASE' 
(Fioc & Rocca-Volmerange 1997). In this paper, we con- 
sider an elliptical galaxy as a host galaxy to relate the 
formation of a QSO. Also, we settle the model param- 
eters so that the color-magnitude relation of present-day 
elliptical galaxies can be reproduced. In this model, we as- 
sume a initial mass function (IMF) as (p = A(m*/MQ)~°-^^ 
for a mass range of [0.1, 60A'/q], where to* is the stel- 
lar mass. The star formation rate (SFR) per unit mass 
, C{t), is assumed to be proportional to the fractional gas 
mass fg{t) = Mg(i)/Mgo, where Mgo is the initial gas mass 



which is lO^^M© and Mg{t) is the total gas mass at time 
t. With incorporated a galactic wind model, the star for- 
mation rate is given by 

C{t) = kfg, {0<t< tw) 

= 0, (i>M, (4) 

where a constant rate coefficient is set to fc=8.6Gyr^^. 
Here, we assume ty, = 7 x 10* yr from the fiducial wind 
model by Arimoto & Yoshii (1987). With this chemical 
evolution model, we can pursue the evolution of the phys- 
ical properties of host galaxy, such as stellar mass, lumi- 
nosity, color, and metallicity. 

3. QSO BH - HOST RELATION 

Based on the present cocvolution model, the evolution 
of the mass of stellar component in the bulge (Mbuige); the 
mass of MDO (Mmdo), and the mass of the supermassive 
BH (-Mbh) are shown in Figure 2, assuming the constant 
Eddington ratio {v = 1). The mass accretion proportional 
to the bulge luminosity leads to the growth of an MDO up 
to 1O*M0, which is likely to form a massive dusty disk in 
the nucleus. However, the matter in the MDO does not 
promptly fall into the BH, because the BH accretion is 
limited by equation(3). The BH mass reaches Mmdo at a 
time tcross- As seen in Figure 2, during t < tcross 

the BH 

mass fraction /bh = MBn/Mbuige increases with time. At 
t > tcross, almost all of the MDO matter has fallen onto 

the central BH, and therefore the BH mass is saturated. 
In the optically thick phase (tthin < t < tvj), the BH frac- 
tion is /bh <C 0.001. In the optically thin phase at t > tw, 
which is the dark gray area {tv, < t < tcross) in Figure 2, 
the BH fraction increases up to /bh — 0.001, which is just 
comparable to the observed ratio. 

It has been argued that color-magnitude relation of 
bulges can be reproduced if a galactic wind sweeps away 
the gas at a wind epoch ty, of several lO^yr (Arimoto & 
Yoshii 1986, 1987; Kodama & Arimoto 1997; Mori et al. 
1997). The evolution of the bulge luminosity is shown in 
Figure 3 with the galactic wind model. Even after the 
galactic wind {t > t^,), Mbh continues to grow until tcross 
and therefore the AGN brightens with time if the Edding- 
ton ratio is constant. After the AGN luminosity (Lagn) 
exhibits a peak at tcross, it fades out abruptly due to ex- 
hausting the fuel of the MDO. As seen in Figure 3, it is 
found that the era of <t< taoss can be divided into two 
phases with a transition time tcrit when ibuigc = -^agn; 
the earlier phase is the host luminosity-dominant phase 
(the dark gray area), and the later phase is the AGN 
luminosity-dominant phase (the light gray area). The life- 
times of both phases are comparable to each other, which 
is about 10*^yr. The AGN-dominant phase is likely to cor- 
respond to ordinary QSOs, but the host-dominant phase is 
obviously different from observed QSOs so far. We define 
this phase as "a proto-QSO". In this phase, /bh rapidly 
increases from 10~^'^ to 10~^'^ in « 10^ years. Also,the 
central massive BH is surrounded by a massive dusty disk 
(> 10^ Mq), which may obscure the nucleus in the edge-on 
view and make a type 2 nucleus. Objects corresponding to 
such host luminosity-dominant proto-QSOs have not been 
identified obscrvationally yet. Hence, the detection of this 
type of QSOs could be a crucial test for the present pic- 
ture. The later fading nucleus could be a low luminosity 
AGN (LLAGN) (e.g. Kawaguchi & Aoki 2001; Awaki et 
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al. 2001). The proto-QSO phase is preceded by a bright 
and optically thick phase, which may correspond to a ul- 
traluminous infrared galaxy (ULIRG) phase. 
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Fig. 2. — BH growth, assuming Mo = 10'' Mq and = 1.0 
[sec equation (3)]. The ordinate is mass in units of Mq. M^uige 
is the mass of stellar eomponent in bulge. Mmdo is the mass of 
MDO. Mbh is the mass of the supermassive BH. iw is the galac- 
tic wind timescale. (cross is defined so that Mmdo = Mbh- The 
light gray area shows the initial optically thin phase in the timescale 
of ~ lO^yr. In this phase, there can exist the massive dusty disk 
(IG^-^M©). In the optically thick phase (tthin < t < tw), the BH 
fraction /bh = A^MDo/A^buige ^ 0.001. In the optically thin phase, 
which is the dark gray area (tw < i < tcross ) , there can be a massive 
dusty disk {> 10*Mq) around a massive BH, and also the BH frac- 
tion is /bh — 0.001, which is just comparable to the observed ratio. 
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Fig. 3. — AGN and bulge luminosity as a function of time. The 
ordinate is the luminosity in units of L©. tcross is the time when 
.^'bulgc = Lagn- Here, we assume that Lagn is the Eddington lu- 
minosity. The phase at t < tw is a bright and optically thick phase, 
which may correspond to a ultraluminous infrared galaxy (ULIRG) 
phase. After the AGN luminosity (Lagn) exhibits a peak at (cross, 
it fades out abruptly. The later fading imcleus could be a low lumi- 
nosity AGN (LLAGN). The optically-thin, bright AGN phase (gray 
area) can be divided into two phases; one is the host-dominant phase 
(proto-QSO), which is the dark gray area (tw < t < tcrit) and the 
other is the AGN-domiiiaiit phase (QSO), which is the light gray 
area (tcrit < * < tcross). The lifetime of both phases are comparable 
to each other, which is about 10**yr. 

The proto-QSO phase may be distinguishable in terms 
of several observable properties like the broad emission 
line width, the color of the host galaxy, and the mctal- 
licity. The time variations of these quantities are shown 
in Figures 4-6. The broad emission line width, which 
corresponds to the virial velocity, Vblr, is assessed by 



an empirical law for the size of broad line region (BLR), 

1 /2 

''BLR = 15^44 light-days (Kaspi et al. 1997), where L44 
is the luminosity at 0.1 — 1/im in units of 10'*'*crg/s. In 
Seyfert galaxies, this relation holds not only for the slow 
growing phase of AGN BHs. but also for NLSyls (Peter- 
son et al. 2000). Provided that the broad line clouds are 
bound in the potential by central BH (Laor et al. 1997) 
and L(0.1 — 1/im) is close to the Eddington luminosity, 
Ltzm ^ 1.2 X 10'*^(MBH/10*MQ)erg/s, the circular veloc- 



ity of broad line clouds (Vblr) is given by 



Vblr — 



GJVL 



BH 



TBLR 



1/2 



1700 



V 1O8M0 



BH 



1/4 



km/s. (5) 



The estimated width of broad emission line is shown in 
Figure 4. It is found that Vblr in the proto-QSO phase is 
less than ~ 1500km/s. This velocity is considerably small 
compared to ordinary QSOs. 
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Fig. 4. — The broad emission line width, Vblr> as a function 
of time. The ordinate is Vblr in unit of km/s. This line rep- 
resents equation (5). In the proto-QSO phase, Vblr is less than 
~1500km/s. 

Also, the colors of host galaxy at the rest and the 
observed bands with assuming the formation redshift of 
Z[ — 4 are shown in Figures 5 (a) and (b), respectively. 
Here, we take into account the effect of dust extinction 
and K-correction, providing Hq = 70km/s/Mpc and the 
Einstein-de Sitter universe. In Figure 5, the thick lines 
show the cases with dust extinction, whereas the thin 
lines denote the cases without dust extinction. At the 
rest frame, as seen in Figure 5 (a), it is found that the U- 
B color is almost identical between proto-QSO and QSO 
phases, because the massive stars responsible for U-B color 
of host galaxy disappear before ~ lO'' yr. However, the 
B-V color in the proto-QSO phase is about 0.5 magni- 
tude bluer than that in the QSO phase. Moreover, the 
V-iT color in the proto-QSO phase is about 0.6 magnitude 
bluer than that in the QSO phase. At the observed frame, 
as seen in Figure 5 (b), the V-K color of the proto-QSO 
phase is about 0.5 magnitude bluer than that in the QSO 
phase, while the noticeable difference is not found in U-B 
and 5- y colors between two phases. We confirm that the 
results do not change significantly for different redshifts 
of galaxy formation (e.g., Zf=3 and 5). This may be an 
important photometric feature that can discriminate the 
proto-QSO phase from the QSO phase. 
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Fig. 5. — (a) U — B, B — V and V — K colors at the rest frame as 
a function of time. The ordinate is the U — B, B — V and V — K 
colors. The effects of dust extinction and K-correction axe taken 
into account. The thick lines show the case with dust extinction. 
The thin lines denote the case without dust extinction. The B — V 
color in the proto-QSOs is about 0.5 magnitude bluer than in that 
of QSOs. Moreover, the V — K color in the proto-QSO is about 
0.6 magnitude bluer than in that of QSOs. (h)U — B, B — V and 
V — K color at the observed frame as a function of time, assuming 
Ho = 70km/s/Mpc, the galaxy formation redshift of 2f=4, and the 
Einstein-dc Sitter model. The upper abscissa shows the redshift. 
The V — K color in proto-QSOs is about 0.5 magnitude bluer than 
the QSO phase. 

As for the metallicity, Figure 6 shows that no appre- 
ciable difference is found between proto-QSO and QSO 
phases. The average mctaUicity of gas in galactic nu- 
cleus, Zblr (the soUd line), is ~ 8Zq, and the aver- 
age metallicity of stars weighted by the host luminos- 
ity, (the dashed line), is ~ 3Zq. These results are 
consistent with the observed metallicity in QSO BLRs 
and the inferred metallicity in present-day elliptical galax- 
ies (Dressier et al. 1984; Hamann & Ferland 1993; 
Shields & Hamann 1997; Dietrich & Whilhelm-Erkens 
2000; Hamann, Foltz & Chaffee 2001; Becker et al. 2001). 
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Fig. 6. — The mean metallicity < Z > as a function of time. The 
ordinate is the metallicity in unit of Zq. The solid line denotes the 
metallicity of gas in galactic nuclei, Zblr- The dotted line shows 
the metallicity of stars weighted by the host luminosity, Z, . The top 
hatched area denotes the observable metallicity in QSO BLR and 
inferred metallicity in present-day elliptical galaxies (Dressier ct al. 
1984; Hamann & Ferland 1993; Shields & Hamann 1997; Dietrich & 
Whilhelm-Erkens 2000; Hamann, Foltz & Chaffee 2001; Becker et 
al. 2001). 

4. A UNIFIED EVOLUTIONAL SCENARIO FOR AN 
ELLIPTICAL GALAXY NUCLEUS 

Based on the results shown in above sections, we at- 
tempt to provide an evolutionary picture for an ellipti- 
cal galaxy nucleus. The present radiation-hydrodynamical 
model is grounded on two key conditions, that is, the star- 
burst event and the optically thick interstellar medium. 
Obviously, starburst galaxies including ULIRGs satisfy 
these two. Starbursts are observed to be often triggered 
by galaxy mergers or interactions (e.g. Borne et al. 2000), 
and also numerical simulations of galaxy mergers actu- 
ally demonstrated the possibility (e.g., Mihos & Hernquist 
1996). Thus, the present model can be considered in the 
context of galaxy mergers. On the other hand, the re- 
cent discovery of high redshift quasars with z > 6 (Fan 
et al. 2001) indicates that the formation of supermassive 
BHs proceeded in appreciably short timescale less than 
lO^yr. In these cases, there may not be enough time for 
well-evolved galaxies to merge to form massive elliptical 
galaxies. However, the present model is still valid if QSO 
hosts satisfy above two conditions. 

Recent high-quality observations suggest that radio 
galaxies, ULIRGs, Lyman break galaxies (LBGs), and Lya 
emitters could bo precursors of spheroidal galaxies at high 
redshifts. As for the radio galaxies, they are categorized as 
radio-loud AGNs. Radio galaxies are thought to have mas- 
sive BHs 10^ Mq), which is comparable to that of the 
radio-loud QSOs, but the ratio of the AGN luminosity to 
the host luminosity {Lagn / Lhuigc) is less than unity (Dun- 
lop et al. 2001). In addition, by the analysis of the spectral 
energy distribution (SED) of high redshift {z > 2) radio 
galaxies, the galactic age is estimated to range from ^ 0.1 
to 2Gyr (Mazzei & Zotti 1996; Pentericci et al. 2001). 
Furthermore, a large amount of dust is detected (e.g., Pa- 
padopoulos et al. 2000), and the gas extended over several 
tens of kpc has solar or supersolar metallicity (Vernet et al. 
2001, Villar-Martin et al. 2001). It is intriguing that all 
these properties of radio galaxies are quite similar to the 
predicted properties of proto-QSOs. In that sense, radio 
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galaxies could be a key candidate for proto-QSOs. 

As for ULIRGs, the X-ray emission (Brandt et al. 1997) 
or Paa lines (Velleux, Sanders & Kim 1999) intrinsic for 
active nuclei have been detected in roughly one forth of 
ULIRGs. Also, recent X-ray observations have revealed 
that ULIRGs with AGN activity mostly have the lower ra- 
tio of the hard X-ray luminosity to the bolometric luminos- 
ity, which is L^/Lboi <C 0.01 (Imanishi et al. 1999; Braito 
et al. 2002). If the AGN luminosity is controlled by the 
Eddington limit, these results indicate that the BH mass 
is considerably small in a ULIRG phase. In the present 
model, if the optically-thick phase (tthin < t < fw) is re- 
garded as a ULIRG phase, the mass ratio MsH/Afbuigo is 
predicted to be much less than 0.001 in the ULIRG phase. 
Also, it is expected that AfBH/A'^buigc grows with metal- 
licity at the later phases of ULIRGs. The present model 
may be a physical picture of the evolution of ULIRGs to 
QSOs proposed by Sanders et al. (1988) and Norman & 
Scovill (1988). 

The present picture also predicts the existence of the 
optically thin phase {t < fthin) before the ULIRG phase 
as shown in Figure 1. This phase has several characteris- 
tic properties; i) The hard X-ray luminosity is relatively 
low Lx = 5 X IO^Lq if I/x = 0. 11/ AGN- ii) The metallicity 
of gaseous component is subsolar, Zblr < 0.1 Zq (Figure 
6). iii) The massive dusty disk with 10^~'^Mq can sur- 
round the galactic center (Figure 2). iv) The seed BH 
with ^ IO'^A/q can form due to the collapse of a rotating 
supermassive star (SMS). Then, the gravitational wave is 
expected to be emitted, which may be detectable by the 
Laser Interferometer Space Antenna (LISA) (de Araujo et 
al. 2001; Saijo et al. 2002). As for the host galaxy, re- 
cently observed high redshift LBGs or Lya emitters may 
correspond to the beginning phase of a bulge (Friaca & 
Terlevich 1999; Matteucci & Pipino 2002). LBGs have 
the luminosity and observed to be optically 

thin and have the metallicity of 0.1 — IZq. Also, they ex- 
hibit strong clustering at z ~ 3. In addition, the chandra 
X-ray observatory has detected the hard X-ray of LBGs 
with the luminosity of ~ 10^ Lq, although it is still uncer- 
tain whether the X-ray emission arises from AGN or not 
(Brandt et al. 2001). Some LBGs have a light profile fol- 
lowing a r^/^ law over a large radial range (Giavalisco et al. 
1996). Moreover, it has been proposed that Lya emitters 
may be the first 10^ yr of the galaxy formation, and there- 
after their luminosity fades by dust rapidly (Malhotra et 
al. 2002). Therefore, the precursor of ULIRGs predicted in 
the present model may correspond to the assembly phase 
of LBGs or Lya emitters. 



5. CONCLUSIONS 

Based on the radiation drag model for the BH growth, 
incorporated with the chemical evolution of the early-type 
host galaxy, we have built up the coevolution model for a 
QSO BH and the host galaxy. As a consequence, we have 
shown the possibility of the proto-QSO phase, which is 
optically-thin and host luminosity-dominant, and has the 
life-time comparable to the QSO phase timescale of a few 
10^ years. We have predicted the observable properties of 
proto-QSOs as follows: (1) The width of broad emission 
line is narrower, which is less than 1500km/s. (2) The BH- 
to-bulge mass ratio, MBn/A^buige, rapidly increases from 
10-5-3 to 10-3-^ in « 10* years. (3) The colors of {B - V) 
at rest bands and {V — K) at observed bands are about 0.5 
magnitude bluer than those of QSOs. (4) In both proto- 
QSO and QSO phases, the metallicity of gas in galactic 
nuclei is Zblr — 8Zq, and that of stars weighted by the 
host luminosity is Z^, ~ 3Zq, which are consistent with 
the observations for QSOs and the elliptical galaxies. (5) 
A massive dusty disk (> lO^M©) surrounds a massive BH, 
and it may obscure the nucleus in the edge-on view to form 
a type 2 nucleus. The predicted properties of proto-QSOs 
are similar to those of radio galaxies. 

The proto-QSO phase is preceded by an optically-thick 
phase before the galactic wind, which may correspond to 
ULIRGs. The present model predicts a low luminosity 
ratio of Lagn to Lboi, which is consistent with the ob- 
served ratio Lx/-^boi ^ 0.01 for ULIRGs. In addition, 
MBn/Afbuige is anticipated to be much less than 10"'^ and 
AfBH/A^buigo grows with metallicity in the ULIRG phase. 

Finally, we can predict the precursor of ULIRGs, which 
is optically thin and their lifetime is ~ 10'' years. This may 
correspond to the assembly phase of LBGs or Lya emit- 
ters. In this phase, the massive dusty disk of « lO^^^A'/0 
exists, the metallicity is subsolar (Z* < O.IZq), and the 
hard X-ray luminosity is Lx ~ 5 x lO^L© if = O.ILagn- 
In addition, the formation of a seed black hole (BH) 
(~ IO^Mq) can occur due to the collapse of a rotating 
supermassive star (SMS) in this phase. Thus, the gravita- 
tional wave may be detectable by the Laser Interferometer 
Space Antenna (LISA). 
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